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ABSTRACT 

Desingularization is the problem of finding a left multiple of 
a given Ore operator in which some factor of the leading co- 
efficient of the original operator is removed. An order-degree 
curve for a given Ore operator is a curve in the (r, d)-plane 
such that for all points (r, d) above this curve, there exists 
a left multiple of order r and degree d of the given opera- 
tor. We give a new proof of a desingularization result by 
Abramov and van Hoeij for the shift case, and show how 
desingularization implies order-degree curves which are ex- 
tremely accurate in examples. 

Categories and Subject Descriptors 

1.1.2 [Computing Methodologies]: Symbolic and Alge- 
braic Manipulation — Algorithms 

General Terms 

Algorithms 

Keywords 

Ore Operators, Singular Points 

1. INTRODUCTION 

We consider linear operators of the form 
L = £o + £id+---+£rd'', 

where £o, ■ ■ ■ ,£r are polynomials or rational functions in x, 
and d denotes, for instance, the derivation ^ or the shift 
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operator x i-^ x + 1. (Formal definitions are given later.) 
Operators act in a natural way on functions. They are used 
in computer algebra to represent the functions / which they 
annihilate, i.e., L • / = 0. 

Multiplication of operators is defined in such a way that 
the product of two operators acts on a function like the two 
operators one after the other: {PL)-f — P{L-f). Therefore, 
if L is an annihilating operator for some function /, and if 
P is any other operator, then PL is also an annihilating 
operator for /. 

We are interested in turning a given operator L into a 
"nicer" one by multiplying it from the left by a suitable P, 
for two different flavors of "nice". First, we consider the 
problem of removing factors from the leading coefficient £r 
of L. This is known as desingularization and it is needed for 
computing the values of / at the roots of £,. (provided it is 
defined there) . Desingularization of differential operators is 
classical [9], and for difference operators, Abramov and van 
Hoeij [2, 1] give an algorithm for doing it. We give below a 
new proof of (a slightly generalized version of) their results. 

Secondly, we consider the problem of producing left mul- 
tiples with polynomial coefficients of low degree. Unlike the 
situation for commutative polynomials, a left multiple PL 
of L may have polynomial coefficients even if P has ratio- 
nal function coefficients with nontrivial denominators and 
the polynomial coefficients of L have no common factors. In 
such situations, it may happen that the degrees of the poly- 
nomial coefficients in PL are strictly less than those in L. 
This phenomenon can be exploited in the design of fast al- 
gorithms because a small increase of the order can allow for 
a large decrease in degree and therefore yield a smaller total 
size of the operator ("trading order for degree"). Degree es- 
timates supporting this technique have been recently given 
for a number of different computational problems [4, 7, 6, 3]. 
Although limited to special situations, these estimates can 
overshoot by quite a lot. Below we derive a general estimate 
for the relation between orders and degrees of left multiples 
of a given operator L from the results about desingulariza- 
tion. This estimate is independent of the context from which 
the operator L arose, and it is fairly accurate in examples. 

2. OVERVIEW 

Before discussing the general case, let us illustrate the 
concepts of desingularization and trading order for degree 



on a concrete example. Consider the differential operator 



we have (go + qid + q2d^)P = ^Qi + Q2, where 



L = -(45 + 25s - 352;^ 



+ 2(33 -9x- 3x' 
(1 + a:) (23 - 20x - x"^ + 2x^)d^ £ <Q[x][d], 



where d 
p := 23 



: That L is desingularizable at (a root of) 
20x - + 2x^ means that there is some other 
operator P G (Q(a:;)[9] such that PL has coefficients in (Q[i] 
and its leading coefHcient no longer contains p as factor. 
Such a P is called a desingularizmg operator for L at p 
and PL the corresponding desingularized operator. In our 
example, 



299 „ 1035 - 104a; - 136a;^ 



P= —d + 
P 



e ^ix)[d] 



is a desingularizing operator for L at p, the desingularized 
operator is 

PL = (-2350 - 2055a; + 104a;^ + 136a;^) 
+ (2151 + 281x + 136a;^)a 

+ (1932 + 931a; - 240a;^ - 136a;^)a^ + 299(1 + x)d^. 

A desingularizing operator need not exist. For example, it 
is impossible to remove the factor a; + 1 from the leading 
coefficient of L by means of desingularization. In Section 3 
we explain how to check for a given operator L and a factor p 
of its leading coefficient whether a desingularizing operator 
exists, and if so, how to compute it. 

Desingularization causes a degree drop in the leading co- 
efficient but may affect the other coefficients of the operator 
in an arbitrary fashion. However, a desingularizing operator 
can be turned into an operator which lowers the degrees of 
all the coefficients. To this end, multiply P from the left by 
some polynomial g e (Q[a;] for which the coefficients of pqP 
have low degree modulo p, i.e., for which qP = -Pi + P2 
where P\,P2 £ (Q[a;][9] and Pi has low degree coefficients. 
In our example, a good choice is g = (—43 + 34a;) /299, i.e. 

Pi = (-22a; + 29) + (-43 + 34a;)9, P2 = 

Since PL has polynomial coefficients, so does 



-PiL = gPL - P2L 
P 

= (-10 - 1652; + 22a;^) + (201 + 65a; - 34a;^)9 
+ (-100 + 109x - 22a:^)a^ - (1 + a;)(43 - 34x)a^ 

This operator has degree deg^(L) + deg^(Pi) — deg^(p) — 
2, compared to deg^(L) + deg^(P) = 3 achieved with the 
original desingularizing operator. There is no left multiple 
of L of order 3 and degree less than 2. There is also none 
of order 4, but there does exist a left multiple of degree 1 
and order 5. It can be obtained from P by multiplying from 
the left by an operator go + gi9 + g2 9^ G (Q[a;][9] of order 2 
for which the coefficients of p^(go + qid + q2d^)P have low 
degrees modulo p"^ : Taking 



Qi = (841 + 580a; - 436x^ - 148a;^ + 59x'' + 122;^ - 'ix'^ 
+ (1697 - 528a; - 752a;^ + 120a;^ + 1272;* - 12a;^ 



Q2 



- 42;'')a + (a; - 7)(-9 + a; + 2a;^)pr + p^d^ 

16(779+3741^) _ 272(69+343:) p. 



89401 



89401 



_ 633+643:-883:^ 
yO 89401 



gi = 



(17a:^+133:-92) 
89401 ' 



'?2 299 



Set Q ■- -^Qi. Then, since PL has polynomial coefficients, 
so does 

QL = (go + qid + q2d^)PL - Q2L 

= (2 + a;) + (-3 + x)d - (8 + 2a;)9^ 
+ (2 - 2x)a^ + (6 + x)d* + (1 + a;)a^ 

Its degree is deg^(L) + deg^((3i) - 3deg^(p) = 1. 

As the factor a; + 1 cannot be removed from L, we cannot 
hope to reduce the degree even further. We have thus found 
that the region of all points (r, d) € such that there is 
a left (Q(2;)[9]-multiple of L of order r and with polynomial 
coefficients of degree at most d is given by ((2,4) + N^) U 
((3,2)+n2)U((5, 1)+n2). 

In Section 4 we explain the construction of the operators 
Q that turn a desingularizing operator into one that lowers 
all the degrees as far as possible, and we give a formula that 
describes the points (r, d) for which such a Q exists. 

3. PARTIAL DESINGULARIZATION 

In this section we discuss under which circumstances an 
operator L admits a left multiple PL in which a factor of 
the leading coefficient of L is removed. This is of interest 
in its own right, and will also serve as the starting point 
for the construction described in the following section. In 
view of this latter application, we cover here a slightly gen- 
eralized variant of desingularization, which not only applies 
to the case where a factor can be completely removed, but 
also cases where only the multiplicity of the factor can be 
lowered. 

Example 1. In the shift case (i.e., dx = [x-^-Vjd), consider 
the operator 

L = (3 + 2;)(9 + 7x + x^) - (33 + 70x -f 47a;^ + 12x=' + x'^)d 
+ {2 + xf{'i + ^x + x^)d'^. 

The factor (2; -|- 2)^ in the leading coefficient cannot he re- 
moved completely. Yet we can find a multiple in which x + 2 
appears in the leading coefficient (in shifted form) with mul- 
tiplicity one only. One such left multiple of L is 

(402 + 208a; + 25x^) - (514 + 743x + 2582;^ + 25a;^)a 
-I- (233 + 378a; + 1832;^ + 25x^)9^ - 9(3 + x)d^ . 

We speak in this case of a partial desingularization. The 
general definition is as follows. We formulate it for operators 
in an arbitrary Ore algebra (D := A[d\ ~ A[d; a, 5] where A 
is a K-algebra (in our case typically A = K[x] or A = K(x)), 
K is a field, a: A ^ yl is an automorphism and S: A ^ A a. 
(T-derivation, i.e., a K-linear map satisfying the skew Leibniz 
rule S{pq) = S{p)q -\- a{p)5{q) for p,q £ A. For any f £ A, 
the multiplication rule in A[9;(t, 5] is df — o-{f)d -\- S{f). 
We write degg(L) for the order of L £ A[d], and if A = 
K[x], we write deg^{L) for the maximum degree among the 



polynomial coefficients of L. For general information about 
Ore algebras, see [5]. 

Definition 2. Let L £ K[a;][9; a, S] and letp G ]K[x] be such 
thatp I lca(L) € K[a:]. We say thatp is removable /rom L at 
order n if there exists some P £ ^{x)\d] with Aegg{P) = n 
and some it), n G K[a::] with gcd(p, w) = 1 such that PL £ 
'K\x][d] and (j-"(lcs(PL)) = :;^^lcs(L). We then call P 
a p-removing operator for L, and PL the corresponding p- 
removed operator, p is simply called removable from L if it 
is removable at order n for some n £ N. 

7/ gcd(p, lca(L)/p) = 1, we say desingulariz [able jingled] 
instead of remov[able\ing\ed], respectively. 

The backwards shift a~" in the definition above is intro- 
duced in order to compensate the effect of the term 9" in P 
on the leading coefficient on L (i.e., lcd{d"L) = cr"(lca(I/)).) 
Moreover, observe that in this definition, removing a polyno- 
mial p does not necessarily mean that the p-removed opera- 
tor has no roots of (some shift of) p in its leading coefficient. 
If L contains some factors of higher multiplicity, as in the 
example above, then removal of a polynomial is defined so as 
to respect multiplicities. Also observe that in the definition 
we allow that some new factors w are introduced when p 
is removed. This is only a matter of convenience. We will 
see below that we may always assume v = w = 1, i.e., if 
something can be removed at the cost of introducing new 
factors into the leading coefficient, then it can also be re- 
moved without introducing new factors. The justification 
rests on the following lemma. 

Lemma 3. Let L £ Vi[x][d; a, 5], let p £ K[x] with p \ 
lcd{L) be removable from L, and let P G K(a;)[9; a, 5] be a 
p-removing operator for L with degg(P) = n. 

l.IfUe lK[x][d] with gcd(lc9(f/),cr"+'*"'5a(^)(p)) = 1, 
then UP ts also a p-removing operator for L. 

2. IfP = Pi -I-P2 for some Pi £ K(x)[9] with degg(Pi) = 
n and P2 £ 'K.[x\\9\, then Pi is also a p-removing op- 
erator for L. 

3. There exists ap-removing operator P' with degg{P') — 
n and with pa~'^ {\cg{P L)) = \cq{L). 

Proof. Let v,w £ K[a;] be as in Definition 2, i.e., gcd(p, w) = 
1 and vpa-'^ilcaiPL)) = wlce(L). 

1. Since PL is an operator with polynomial coefficients, 
so is UPL. Furthermore, with u — leg (17) and m — degg{U) 
we have 

vpa~"-"^{\ca{UPL)) = o--"-™(«)wlce(L). 

Since gcd{u,a"+"'{p)) = 1, we have gcd(cr-"-™(u)w,p) = 
1, as required. 

2. Clearly, P2 £ Vi[x][d] implies P2L £ K[x][d]. Since also 
PL £ K[x][d], it follows that 

PiL = (P - P2)L = PL~ P2L £ ^[x\[d]. 

If degg(P2) < n, then we have \cg{PL) — \cq{PiL), so there 
is nothing else to show. If degg(P2) = n, then leg (Pi L) = 
Ice (Pi) ^ lca(P2L) and therefore 

vpa-"{lca{PiL)) = i;po--"(lcs(PL) - lc9(P2L)) 
= (w - vpa-"{\ce{P2)))lce{L). 



Since gcd{p,w ~ vpa~'^{\cg{P2))) = gcd(p, w) — 1, the claim 
follows. 

3. By the extended Euclidean algorithm we can find s,t £ 
K[a;] with 1 — sw -\- tpv. Then a"{s)P is p-removing of 
order n by part 1 (cr"(s) is obviously coprime to ct"(p)), and 
its leading coefficient is 

^„/£U,N 1 
\ pv / a"(pv) 

By part 2 we may discard the polynomial part cr'^{t), ob- 
taining a p-removing operator P' with the desired property. 
■ 

The lemma implies that if there is a p-removing operator 
at all, then there is also one in which all the denominators 
are powers of f"(p) (because any factors coprime with p 
can be cleared according to part I), and where all numera- 
tors have smaller degree than the corresponding denomina- 
tors (because polynomial parts can be removed according to 
part 2). 

Similarly as in the proof of part 3, we can also reduce the 
problem of removing a composite polynomial to the problem 
of removing powers of irreducible polynomials. For exam- 
ple, if p = P1P2 is removable from L, where pi,P2 £ lK[x] 
are coprime, then obviously both pi and P2 are removable. 
Conversely, if pi and P2 are removable, and if Pi, P2 are re- 
moving operators of orders ni, n2 with lc8(Pi) = l/a"^ (pi) 
and lca(P2) = 1/(t"^(p2), then for n = max{ni,n2} and 
ui,U2 £ K[a;] with 

UiCt"{p2) + U2a"{pi) = gcd((T"(pi),(j"(p2)) = 1 

the operator P := md"'"^ Pi 4- U29""""P2 G lK(a;)[9] is 
such that PL £ K[x][d] and Ice (PL) = lce(L)/cr"(p). 

In summary, in order to determine whether a polynomial 
p = Pi^P2^ ■ ■ -p^ is removable from an operator L, it suf- 
fices to be able to check for an irreducible polynomial pi 
and a given ki > 1 whether p*'' is removable. Let now p 
be an irreducible polynomial and fc > 1. If there exists a 
p'^-removing operator, then it can be assumed to be of the 
form 

p = P° + d+- ■ .+ ^"-^ ^ 9" 

cr"(p)<=r) cr"(p)'^i (T"(p)'="-i a"{p)'' ' 

for some eo,...,e„_i G N, and po, . . . , Pn-i G K[a;] with 
deg^(pi) < eideg^(p). In order to decide whether such an 
operator exists, it is now enough to know a bound on n as 
well as a bound e on the exponents e,;, for if n and e are 
known, we can make an ansatz pi = TTjZ^Vi.j^^ with un- 
determined coefficients Pi.j, then calculate PL and rewrite 
all its coefficients in the form a/a"{pY -\-b for some polyno- 
mials a, b depending linearly on the undetermined Pi j , then 
compare the coefficients of the various a's with respect to x 
to zero and solve the resulting linearly system for the pij. 

How the bounds on n and e are derived depends on the 
particular Ore algebra at hand. In this paper, we give a 
complete treatment of the shift case {{ap){x) = p{x 4- 1), 
5 = 0) and make some remarks about the differential case 
(ct = id, 5 = ^). For other cases, see the preprint [8]. 

3.1 Shift Case 

In this section, let K[x][9] denote the Ore algebra of re- 
currence operators, i.e., a is the automorphism mapping x 
to a; -|- 1 and S is the zero map. This case was studied by 



Abramov and van Hoeij [2, 1]. We give below a new proof 
of their result, and extend it to the case of partial desingu- 
larization. For consistency with the differential case, we for- 
mulate the result for the leading coefficients, while Abramov 
and van Hoeij consider the analogous for the trailing coeffi- 
cients. Of course, this difference is immaterial. 

We proceed in two steps. First we give a bound on the 
order of a removing operator (Lemma 4), and then, in a 
second step, we provide a bound on the exponents in the 
denominators (Theorem 5). As explained above, it is suffi- 
cient to consider the case of removing powers of irreducible 
polynomials, and we restrict to this case. 

Lemma 4. Let L = h) + hd + h Ird'' G K.[x\[d] with 

£o,£r 7^ 0, and let p be an irreductble factor oflcg{L) such 
that p'' is removable from L for some k > 1. Let n G N 6e 
s.t. gcd (cr" (p), £()) 7^ f and gcd{a"^ (p) , £o) = f for all m> n. 
Then p^ is removable at order n from L. 

Proof. By assumption on L, there exists a p*-removing op- 
erator P, say of order m, and by the observations following 
Lemma 3 we may assume that 

PO , Pi Q , , Pm 



P = 



(7'"(p)=0 (J™(p)'=l 



-9" 



for ei G N and pi G K[a::] with deg^(pi) < eideg^(p) (z = 
0, ...,m). We may further assume gcd(o-'"(p),pi) = f for 
i = Q, . . . ,m (viz. that the Ci are chosen minimally). 

Suppose that m > n. We show by induction that then 
eo = ei = • • ■ — Cm-n-i = 0, so that pi = for i = 
0, . . . , m — 71 — f , i.e., the operator P has in fact the form 

Pm — n 



P = 



-d" 



cr'"(p)'= 



-a" 



Thus 9" ""P G K(a::)[9] is a ^''-removing operator of or- 
der n. 

Consider the operator T ■- YJlt^ ^id" — PL G K[x][d]. 
From 



Hp) 



—lo ~ t{) G K[a;] it follows that eo = 0, because 
jcd(a'"(p),po) = gcd(a'"(p),£o) = 1 



by the choice of po and the assumption in the lemma, re- 
spectively, and this leaves no possibility for cancellation. 

Assume now, as induction hypothesis, that eo = ei = 
• ■ ■ = Ci-i — for some i < m — n. Then from 



Pi ill, N , P»-l 



a^{pY 

Pi 
a^{pY 



a^{py 



+ 



PO 



it follows that a'^{pY' \ pia^{£o)- By the choice of pi we have 
gcd(a'"(p),pi) = f and by the assumption in the lemma 
we have gcd((T™'~*(p), ^o) ~ 1 (because m ~ i > n), so it 
follows that Ci = 0. Inductively, we obtain eo = ei = ■ • ■ = 
Cm-n-i ~ 0, which completes the proof. ■ 

It can be shown that p cannot be removed from L if a" (p) 
is coprime with the trailing coefficient of L for all n G N by 
a variant of [2, Lemma 3.], so the above lemma covers all 
situations where removing of a factor is possible. 

In order to formulate the result about the possible ex- 
ponents in the denominator, it is convenient to first intro- 
duce some notation. Let us call two irreducible polynomi- 
als p, 5 G K[a;] \ {0} equivalent if there exists n £ W, such 
that a"{p)/q G K. We write [q] for the equivalence class 



of g G IKlx] \ {0}. If p,q are equivalent in this sense, we 
write p < q a a"{p)/q G K for some n > 0, and p > q 
otherwise. 

The irreducible factors of a polynomial u G K[a;] can be 
grouped into equivalence classes, for example 

u = (x — 4){x - lfx{x + if {2x - 5){2x + 3f(2x + 9) 

X [x^ + 5x + l){x^ + llx + 25f . 
■■ « ' 

For any monic irreducible factor p of it G K[a::], let Vp{u) 
denote the multiplicity of p in u, and define 

u<p(u) := m&x{vq{u) | g G [p] : p > g }. 

For example, for the particular u above we have Vx-4{u) — 1, 
v^x-4{u) = 0, v^x+i{u) = 3, and so on. 

Besides being applicable not only to desingularization but 
also removal of any factors, the following theorem also refines 
the corresponding result of Abramov and van Hoeij in so far 
as their version only covers the case of desingularizing L at 
some p with w>p(lca(L)) = whereas we do not need this 
assumption. 

Theorem 5. Let L = £o + £18 -\ h £rd'' G TK[x][d] with 

£0, £r 0, and letp be an irreducible factor of £r such thatp*' 
is removable from L for some k >1. Let n G N 6e such that 
gcd(cr"(p),4) / 1 and gcd(a'"(p), 4) = 1 for all m > n. 
Then there exists a p*" -removing operator P for L and p of 
the form 



Po 



+ ■ 



Pi 



-d+--- + 



Pn 



a"{pY 



dn 



(7"(p)eo cr"(p)'^ 

for some e^ G N and pi G K[a;] with 

1. deg^(pi) < eideg^(p) anrf gcd(cr"(p),pi) = 1, and 

2. Ci < k + nv^p(lca{L)) 

for i = 0, . . . ,n — 1, and pn — 1, e„ — k. 

Proof. Lemmas 3 and 4 imply the existence of an operator P 
with all the required properties except possibly the exponent 
estimate in item 2. Let P be such an operator, and consider 
the operator T — Yli^o tid" := PL G K[x][d]. 

Let e = max{ei, . . . , e„} and P ~ Yli^oPid'' ■= (^"{pTP- 
Then p, = a"{py-''^pi {i = 0, . . . , n) and a"{pyT = PL 
and gcd(po, . . . ,Pn, cr"{p)) = 1. 

Abbreviating v :— v^p{\cd{L)), assume that e > k + nv. 
We will show by induction that then pi contains o-"(p) with 
multiplicity more than iv for i = n,n — I, ... ,0, which is 
inconsistent with gcd(po, • • • ,Pn, cr"(p)) = f . 

First it is clear that p„ = a" (pYpncr" (£r) contains cr"(p) 
with multiplicity > e — fc > nv, because P is p^-removing. 
Suppose now as induction hypothesis that there is an i > 
such that (t" (p)-* | pj for j — n,n — l, . . . ,i + l. Consider 
the equality 

Cr"(p)''ti + r = Pi(j\tr) + Pi + ia''^^{£r-l) H hPn^r-n, 

where we use the convention £-j ~ Q for j < 0. The induction 
hypothesis implies that (t"(p)''+^'"+^ | pj for j — n,n — 
I,...,i + I. Furthermore, since (i + l)v < nv < e, we 
have cr"(p)('+^'"+^ | a"{pYti+r. Both facts together imply 
cr"(p)(»+i)''+i I pia'{£r). The definition of v ensures that 
(t"(p) is contained in a^{£r) with multiplicity at most v, 
so it must be contained in pi with multiplicity more than 
{i + l)v ~ v — iv, as claimed. ■ 



3.2 Differential Case 

In this section K[x][9] refers to the Ore algebra of differen- 
tial operators, i.e., cr = id and 5 = Let L £ K[a;][9] and 
suppose for simplicity that p = a; is a factor of \cd{L). In [1], 
the authors show that L can be desingularized at x if and 
only if a; = is an apparent singularity, that is, if and only if 
L{y) — admits degg{L) linearly independent formal power 
series solutions. The authors furthermore give an algorithm 
to find an operator P such that if ^ is either an ordinary 
point of L or an apparent singularity of L, then ^ is an or- 
dinary point of PL. Therefore this algorithm desingularizes 
all the points that can be desingularized. The authors also 
give a sharp bound for degg(P). The authors furhtermore 
give some indications concerning partial desingularizations. 
It would be interesting to give a complete algorithm for par- 
tial desingularizations. 

4. ORDER-DEGREE CURVES 

We now turn to the construction of left multiples of L with 
polynomial coefficients of small degree, and to the question 
of how small these degrees can be made. As already in- 
dicated in Section 2, we start from an operator P which 
removes some factor from the leading coefficient of L, say it 
removes a polynomial p of degree k. According to Lemma .3, 
we may assume that lca(P) = l/a'^''^^^^\p) and that all 
other coefficients of P are rational functions whose numera- 
tors have lower degree than the corresponding denominators. 
Thus we already have deg^(PL) < deg^(L) — 1. Further- 
more, if q is any polynomial with deg^(q) < deg^(p) = k, 
then multiplying P by q (from left) and removing polyno- 
mial parts by Lemma 3.2 gives another operator Q with 
deg3;(Qi) < deg^(L) — 1. All the operators Q obtained in 
this way form a K-vector space of dimension k. Within this 
vector space we search for elements where deg^{QL) is as 
small as possible. Forcing the coefficients of the highest de- 
grees to zero gives a certain number of linear constraints 
which can be balanced with the number of degrees of free- 
dom offered by the coefficients of q, as illustrated in the 
figure below. As long as we force fewer than k terms to 
zero, we will find a nontrivial solution. 



If we want to eliminate k terms or more in order to get 
a result of even lower degree, we need more variables. We 
can create k more variables if instead of an ansatz qP we 
make an ansatz (go -I- qid)P for some qo,?i £ TK[x] with 
deg^(go), deg^(gi) < k. Again removing all polynomial parts 
from the rational function coefficients we obtain a vector 
space of operators Q with deg^{QL) < deg^(L) — 1 whose 
dimension is 2k. The additional degrees of freedom can be 
used to eliminate more high degree terms, the result be- 
ing an operator of lower degree but higher order. If we let 
the order increase further and for each fixed order use all 
the available degrees of freedom to reduce the degrees to 
minimize the degrees of the polynomial coefficients, a hy- 
perbolic relationship between the order and the degree of 
QL emerges. In Theorem 9 below, we make this relation- 
ship precise, taking into account that for a given operator L 



the leading coefficient may contain several factors p that are 
removable at different orders n. The resulting region of all 
points (r, d) £ for which there exists a left multiple of L 
of order r with polynomial coefficients of degree at most d 
is then given by an overlay of a finite number of hyperbolas. 

Before turning to the proof of this theorem, let us illus- 
trate its basic idea with the example operators from Sec- 
tion 2. 

Example 6. Let L G (Q[a;][a], p G (Qlx], and P G (Q(a;)[9] be 
as in Section 2. Recall that p is an irreducible cubic factor of 
\cg[L) and that P is a p -removing operator for L. We have 
P — ^9+ ^ for some pi,po G <Q[x] with deg^(pi) = and 
deg^(po) ~ 2. We have seen in Section 2 that there is an 
operator Q G (Q(a;)[9] of order 3 such that QL G (Q[a:][9] and 
deg^{QL) = 1. Our goal here is to explain why this operator 
exists. 

Make an ansatz Qi = (go + qid + q2d^)P with undeter- 
mined polynomials go , gi , <?2 ■ After expanding the product 
and applying commutation rules, Qi has the form 

Pi 92 gs ^ (. . .)g2 + (■ ■ Qgi g2 

p p2 

, (. . .)g2 + (. . .)gi + (. . .)go ^ , (. . .)g2 + (. . .)gi + (. ■ .)go 

where the (. . .) are certain polynomials whose precise form 
is irrelevant for our purpose. 

Note that by Lemma 3.1, QiL G (Q[2;][9] regardless of the 
choice of go , gi , g2 , and that by Lemma 3. 2, this property is 
not lost if we add to Q\ some operator in (Q[a;][9] of our 
choice. Therefore, if Q2 G (Q[a;][c)] is the operator obtained 
from p'^Qi G (Q[a;][9] by reducing all the coefficients mod- 
ulo p^, then p~'^Q2L G (Q[a;][9], still regardless of the choice 

ofqo,qi,q2. 

The coefficients of Q2 depend linearly on the undeter- 
mined polynomials go,gi,g2. If we choose their degree to be 
deg^(p) — 1 = 2, then we have 3(2 + 1) = 9 variables for the 
coefficients o/go, gi, g2. Choosing a higher degree would give 
more variables but also introduce undesired solutions such as 
go = gi = g2 = p, for which the reduction modulo p^ leads to 
the useless result Q2 = 0. This cannot happen if we enforce 
deg^('70 < deg^(p). 

The operator p ^Q2L has degree 

deg^(Q2) +deg^(L) - 3deg^(p) = deg^(g2) - 5, 
which is equal to 1 if deg^{Q2) — 6. A priori, the degree of 



Q2 in x may be up to deg^(p ) — 1 = 8. In order to 
it down to 6, we equate the coefficients of x^d-' for i = 7,8 
and j — 0, ... ,3 to zero. This gives 8 equations. As there 
are more variables than equations, there must be a nontrivial 
solution. 

For formulating the proof of the general statement, it is 
convenient to work with an alternative formulation of remov- 
ability, which is provided in the following lemma. Through- 
out the section, K[a;][9] = K[a;][9;a, 5] is an arbitrary Ore 
algebra. 

Lemma 7. p G K[a::] is removable from L G K[a::][9] at or- 
der n if and only if there exists P G K[x][9] with degg(P) = 
n and PL G a"{p)lcg{P)K[x][d]. 



Proof. 



Po = 



Hp) ica{-P) 



P is a p-removing operator. 



Start from a p-removing operator of the form 

n-l 

Pi 



i=0 



cr"(p) = 



cr"(p) 



dn 



and set P = a"{pyPo where e — maxjeo, . • • , e„„i, 1} > 1. 
Because of PqL £ 'K[x][d] it follows that 

PL G a"{pyK[x][d] = o-"(p)lcs(P)K[x][9]. ■ 

The next lemma is a generalization of Bezout's relation 
to more than two coprime polynomials, which we will also 
need in the proof. 

Lemma 8. Let ui, . . . , Um £ K[a::] be pairwise coprime and 
u — uiU2---Um, and let vi, . . . ,Vm £ IK[a:] be such that 
deg^(««) < deg^(ui) (i = 1, . . . ,m). If 



Eu 
Vi — 



then vi = V2 



Vm = 0, 



Proof. Since the Ui axe pairwise coprime, Ui \u/ui for all i. 
However, Ui \ u/uj for all j ^ i. Both facts together with 
51^™,! Viu/ui — imply that Ui \ Vi for all i. Since deg^{vi) < 
deg^(iti), the claim follows. ■ 

Theorem 9. Let L £ 'K[x][d], and let pi, . . . ,pm G TK[x] be 
factors oflcg{L) which are removable at orders ni, . . . ,n„i, 
respectively, so that the cr"* (pi) are pairwise coprime. Let 
r > degQ{L) and 



d>degJL) 



degs(L) + l 



where we use the notation (x)^ := max{a::,0}. Then there 
exists an operator Q G ]K(a;)[9]\{0} such that QL G K[a;][9] 
and degg(QL) = r and Aeg^{QL) — d. 



Proof. Let r > degg(L), and set s 



degg{L) so that 



s — degg(Q). We may assume without loss of generality 
that s is such that 1 . r > , , = 1 



r-dcga(L) + l 



^ > for all i 



by simply removing all the Pi for which 1 — < from 
consideration. We thus have s > Ui for all i. 

Lemma 7 yields operators Pi G K[a;][9] of order Ui with 
PiL G a"^{p^) lciP,)K[x][d]. Set 



m s — rii 



{pi)a\li 



where U = \cg(Pi). Consider the ansatz 



Also observe that any operator Q2 obtained in this way is 
nonzero unless all the are zero, because if k is maximal 
such that at least one of the g^^fc is nonzero, then 



lce(Qi) = ^^qi, 



is nonzero by Lemma 8. Furthermore, lca((5i) ^ mod 
q because deg^(gi,fc) < deg^(pi) implies deg^(lca(Qi)) < 
deg^('7)- 

The ansatz for the g^j gives X^i" + cleg^(pi) vari- 

ables. Plug this ansatz into Q\ and reduce all the polynomial 
coefficients modulo g, obtaining an operator Q2 of degree less 
thandeg^(g) = EIli(s-"'i + l)(deg^(pO+deg\r(^0)- Then 
for each of the s + 1 polynomial coefficients in Q2 equate the 
coefficients of the terms x-' for 



> ^(s - nO(deg,(pO +deg,aO) + 



to zero. This gives altogether 



E"i^.deg,(pO 
s + 1 



+ 1) ( E((^-'^> + l)(degJpO+degJ/0) - 1 - deg.aO) 



EI^i'^ideg^(pO 



^(s-ni)(deg^(pi)+deg^(Zi)) ^ 

i = l 

:(s + l)^degjp0-l- 



s + 1 

E^^i^i deg^(pO 
s + 1 



equations. The resulting linear system has a nontrivial so- 
lution because 



#vars — #eqns 

m 

= ^(s -n, -f l)deg^(p,) 



-(s + l)(^degjp0-l- 

m / 

= -^nideg^(pi)-(s + l)(-l 



EIli^» deg^(pi 
s + \ 

EI^i"-» deg^(pO 



s + \ 



> - ^nideg^(pi) + j^X]'^»'^^S:r(P> 



a^+"'(pi)aHli) 



d'Pi 



for undetermined polynomial coefficients qij (i — 1, . . . ,m; 
J = 0, . . . , Wi) of degree less than deg^(pi). Regardless of the 
choice of these coefficients, we will always have Qi G K[a;][9] 
and QiL G gK[x][9]. Also, for arbitrary R G K[a::][9] and 
Q2 = Qi- qR we have Q2 G TK[x][d] and Q2L G qTK[x][d]. 
This means that we can replace the coefficients in Qi by 
their remainders upon division by g without violating any 
of the mentioned properties of Qi. 



By construction, the solution gives rise to an operator Q2 G 
K[a;][9] of order at most n with polynomial coefficients of 
degree at most 



- n-i)(deg^(pO + deg^(/i)) 



T,T=i'ni deg^(pO 
s + 1 



for which Q2L G gK[a;][a]. Thus if we set Q = iQ2 G 
K(a;)[S], we have dega(QI/) = degg{L)+s = r anddeg^(QI/) 



is at most 

degJL) + deg^(Q2) -degjg) 

m 

< deg^iL) + Y,{s - n0(degjp0 + degj/.)) 

i=l 

EIli"-» deg^(Pi 
s + 1 

- ^(s - n. + l)(deg^(pO + deg^(/0) 
1=1 



+ 



< deg JL) - ^ deg JpO + 



s + 1 



^ degJL) 



■ rri. 



as required. (The final step uses the facts [— xj = — [a;] and 
\x + n] = [a;] + n for a; £ E, and n € Z.) ■ 

Example 10. 1. Consider again the example from Sec- 
tion 2. There we started from an operator L £ K[a;][9] 
of order 2 and degree 4 for whtch there exists a desin- 
gularizmg operator P of order 1 which removes a poly- 
nomial p of degree 3. According to the theorem, for 
every r > 2 exists an operator Q £ K[a;][9] with QL £ 
K[a;][9], degg(QL) < r and 



d--deg^{QL) <4 



r-2 + 1/ 



3 = 



This hyperbola precisely predicts the order-degree pairs 
we found in Section 2: 



r 


2 


3 


4 


5 6 


7 8 9 


d 


4 


2 


2 


1 1 


1 1 1 



2. Consider the sequence (a„)55^o defined by 
r(2n + fc)r(n - fc + 2) 

k 



r(2n-fc)r(n + 2fc) 



(n £ N). 



Zeilberger's algorithm finds an annihilating operator L 
of the form 

9(1 + n)(l + 3n)(2 + 3n)^(3n + 4)p(n + 1) 
+ (. . . degree 16. . . )9 + (. . . degree 15. . . 

- 10n(8 + 5n)(9 + 57i)(ll + 5n)(12 + 5n)p(r^)a^ 

where d represents the shift operator and p is a certain 
irreducible polynomial of degree 10. This polynomial 
is removable of order 1. Therefore, by the theorem, we 
expect left multiples of L of order r and degree bounded 
by 

■ - 2 



16 



1 

3 + 1 



10 : 



In the figure below, the curve d — (solid) is con- 

trasted with the estimate d = (dashed) derived 

last year for this example [6] as well as the region of 
all points (r, d) for which a left multiple of L of order r 
and degree d exists (gray). The new curve matches pre- 
cisely the boundary of the gray region, even including 
the very last degree drop (which is not clearly visible 



on the figure): for r 
r — 13 we have 



6r-2 
r-2 



- 12 we have 
■■ 6.9 < 7. 



&r-2 

r-2 



7 and for 



30 




3. Consider the minimal order telescoper L for the hyper- 
exponential term in Example 15.2 in [7]. It has order 3 
and degree 40. The leading coefficient contains an irre- 
ducible polynomial p of degree 23 at order 1 and other- 
wise only non-removable factors. Theorem 9 therefore 
predicts left multiples of L of degree r and degree 

17r - 11 



40 



1 

r-3 + 1 



23 : 



r-2 



for all r £ TM. Again, this estimate is accurate, while 
the estimate derived in [7] overshoots. 



50 r 



30 



10 




4. Operators coming from applications tend to have lead- 
ing coefficients that contain a single irreducible polyno- 
mial of large degree which can be removed at order 1, 
besides factors that are not removable. But Theorem 9 
also covers the more general situation of factors that 
are only removable of higher order, and even the case 
of several polynomials that are removable at several or- 
ders. As an example for this general situation, con- 
sider the operator 

L = 8(l+T)(l+22:)^(37+32)^(14+32a;+26a;^+7x=')^ 

-9(l+3x)''(2+3a;)^(l+x+5x^+7x^)^a, 

where d represents the shift operator. From its lead- 
ing coefficient, the polynomial (1 + a; + bx^ + Tx'^Y is 
removable at order 1, and in addition, (l + 3a;)^ is re- 
movable at order 12. The remaining factors are not 
removable. According to Theorem 9 we expect that L 
admits left multiples of order r and degree 



32 - 21 



for all r £ N. It turns out that this prediction is again 
accurate for every r. Observe that in this example the 
curve is a superposition of two hyperbolas. 



d 




"o 5 10 15 20 25 '' 

In conclusion, we believe that removable factors provide 
a universal explanation for all the order-degree curves that 
have been observed in recent years for various different con- 
texts. We have derived a formula for the boundary of the 
gray region associated to a fixed operator L, which, although 
formally only a bound, happens to be exact in all the ex- 
amples we considered. This does not immediately imply 
better complexity estimates or faster variants of algorithms 
exploiting the phenomenon of order-degree curves, because 
usually L is not known in advance but rather the desired 
output of a calculation, and therefore we usually have no in- 
formation about the removable factors of \cq{L). However, 
we now know what we have to look at: in order to improve 
algorithms based on trading order for degree, we need to 
develop a theory which provides a priori information about 
the removable factors of lca(L). In other words, our re- 
sult reduces the task of better understanding order-degree 
curves to the task of better understanding what causes the 
appearance of removable factors in operators coming from 
applications. 
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